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factors, and their interpretation in terms of missing pion-cloud effects. 
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1. Introduction 

Probing hadrons with electromagnetic, axial and pseudoscalar currents reveals their basic 
structure properties and provides a connection with the underlying quark and gluon dynamics in 
Quantum Chromodynamics (QCD). While the nucleon's axial structure is experimentally more dif- 
ficult to access, an abundance of information has been collected for photon-induced processes that 
are described by A^A'^y elastic and A^Ay transition form factors. Precision measurements have stim- 
ulated the development of tools to address questions related to quark orbital angular-momentum 
correlations in the perturbative domain, the transition between perturbative and non-perturbative 
regions, or pion-cloud rescattering effects in the chiral and low-momentum region. The associated 
chiral non-analyticities stemming from the nucleon's 'pion cloud' have been frequently discussed 
when connecting results from lattice QCD, chiral effective field theories and quark models with 
experiment. 

A complementary framework for studying hadron phenomenology is the one via Dyson- 
Schwinger equations (DSEs). They interrelate QCD's Green functions and provide access to non- 
perturbative phenomena such as dynamical chiral symmetry breaking and confinement, see ^ 
for reviews. The investigation of hadron structure in the Dyson-Schwinger approach proceeds via 
covariant bound-state equations, i.e., the Bethe-Salpeter equation (BSE) for mesons and the covari- 
ant Faddeev equation for baryons [^, Q]. The approach has several benefits: it is Poincare-covariant 
throughout every step and provides access to all momentum scales and all quark masses without the 
need for extrapolations. Since one operates directly with QCD's degrees of freedom, observable 
phenomena at the hadron level can be systematically traced back to their microscopic origin. 

The drawback of the approach is its necessity of truncations. Owing to the numerical complex- 
ity of the Faddeev equation, present baryon calculations have been performed in a rainbow-ladder 
truncation, where qqq interactions are neglected and the qq and qq interactions are modeled by 
a dressed gluon exchange. As a consequence, several phenomenologically important features are 
missed in the resulting form factors. A characteristic example is the absence of pion-cloud contri- 
butions in their chiral and low-momentum structure. The relevant gluon topologies that generate 
pion-cloud effects at the hadron level are not captured by a rainbow-ladder truncation which there- 
fore represents the baryon's 'quark core'. In the case of the A^A/ transition form factors discussed 
below, an additional quark-diquark simplification is made, where scalar and axialvector diquark 
correlations approximate the qq scattering matrix and lead to an effective two-body description. 

In the following we will summarize recent results for the nucleon's electromagnetic, axial 
and NAj transition form factors. More detailed discussions, result tables as well as references to 
experimental and lattice data which are frequently used in the plots for comparison can be found in 
Refs. [1-0]. 

2. Covariant Faddeev approach 

The description of baryon structure properties in the Dyson-Schwinger approach requires 
knowledge of the nucleon and A bound-state wave functions and their microscopic ingredients 
in terms of QCD's Green functions. A convenient starting point is given by the three-quark con- 
nected and amputated scattering matrix T. It encodes the relevant information on baryons which 
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Figure 1: Covariant three-body equation for a baryon wave function, of. Eq. (2.3). 

correspond to poles in T. At a given pole for a baryon with mass M the scattering matrix assumes 
the form _ 

T ~ ^ > ^ ^ (2.1) 

which defines the baryon's covariant wave function *F, and *F is its charge conjugate^ . 

Of course, if the scattering matrix were known a priori, the masses of nucleon and A together 
with their wave functions could be directly extracted. This is not the case and thus we aim for rela- 
tions that allow to circumvent the explicit determination of T in practical applications. T satisfies 
a scattering equation, i.e., the nonperturbatively resummed Dyson series 

3 

T = K + KGoT, with Go = S(8)S(g)5 and K = %] + £ ^^^^j ® . (2.2) 

By construction, the kernel K is the sum of a three-quark irreducible contribution A'[3] and per- 
muted two-quark irreducible kernels ^(^j, where the subscript a stands for the respective spectator 
quark. S denotes the dressed quark propagator. The combination of Eqs. (^?T|-p!2[), evaluated at a 
bound-state pole = — M^, yields a self-consistent integral equation for the baryon wave function, 
cf.Fig.|T} 

^' = KGq^. (2.3) 

It can be solved once the dressed quark propagator and the qq and qqq kernels, which encode the 
interactions at the quark-gluon level, are determined. Naturally, all these relations are equally valid 
in the meson case if the three-quark scattering matrix and kernel are replaced by their qq analogues 
and Go is taken as the qq propagator product. 

The coupling of the baryon to an external qq current, on the other hand, is reflected by the 
'gauged' scattering matrix T'' whose residue at the bound-state pole defines the current matrix 
element : _ 

. — I y I ! (2 4) 

{Pj+M^)(Pf+M^)' 

Depending on the type of current, which we generically denote by the index }JL, the respective 
matrix element contains for example the electromagnetic, axial or pseudoscalar form factors of 
the baryon. Here Pi and Pf are the incoming and outgoing baryon momenta and and *Pj are 
the corresponding wave functions. They need not describe the same type of baryon; for instance. 



Eq. ( 2.4 ) could also be applied to the Nhy transition. 



The requirement that the current couples linearly to all internal building blocks of the scattering 



matrix T implies that it has the formal properties of a derivative. Eq. ( |2.2| ) can then be used to 



'For simplicity, we will use the term 'wave function' here both for the bound-state amplitude ^ and the bound-state 
wave function GqI', where the latter has quark propagator legs attached. 
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Figure 2: General expression for a baryon's current matrix element given in Eq. (2.5). The qq 
vertex, dressed quark propagator, and qq and qqq kernels are sandwiched between incoming and 
outgoing baryon wave functions. 



resolve T'^ to a coupling to the dressed quark propagator and the kernel K. The general expression 
for a baryon's non-perturbative current is thereby obtained as follows [||, ^, ||]: 



(2.5) 



It is illustrated in Fig. ^ and consists of an impulse-approximation diagram and further contributions 
involving the qq and qqq kernels. Its ingredients are given by 

= I rf,) ® , = £ rf,, ^k^.^ + £ ^^-^ ^k^^+k^,,, (2.6) 

£1=1 u=\ a=\ 

where the quark labels {a,b,c] are an even permutation of {1,2,3}. 

Let us characterize the external current (for example, electromagnetic, axialvector or pseu- 
doscalar) by Fq e {Z2/7^, 22/57^, Z^iy$ }, equipped with appropriate flavor structures, where Z2 
and Z4 are renormalization constants. The microscopic coupling of the current to the quark is then 
represented by the respective qq vertex F'^ which satisfies a Dyson-Schwinger equation: 



(2.7) 



where T denotes now the qq scattering matrix and K the qq kernel. Pictorially speaking, this 
amounts to the sum of a pointlike part plus all possible reaction mechanisms between quark and 
antiquark which constitute the scattering matrix. In the second step we have exploited the scattering 



equation (2.2) for T to obtain an inhomogeneous Bethe-Salpeter equation for the vertex which, in 



analogy to the bound-state equation ([2.3|), allows to determine the vertex self-consistently from the 
qq kernel K. 

The appearance of the quark-antiquark T-matrix in the defining equation for F'' entails that 
the vertex contains meson poles whenever the bare structure Fq has non- vanishing overlap with the 



respective meson wave function, cf . Eq. (2.1) 



M 



'M 



Tr /»FmGoF; 



(2.8) 



Here, Q = Pf — Pi is the total qq momentum that flows into the vertex and the index 'M' stands 
for meson. Since F'^ enters the form factors diagrams via Eq. ( |2.6| ) it is clear that these poles must 
also appear in the timelrke structure of the form factors where they set the relevant scales: the 
p— meson dominates electromagnetic processes, the axialvector meson ai appears in axial (isovec- 
tor) form factors, and the pion and its excitations in pseudoscalar form factors of hadrons. 
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The same reasoning was generalized in Ref. [10] to derive the hadron's coupling to two exter- 
nal currents with qq quantum numbers: 

P}=P}=-M^ ^Ifp^"^. 



{P}+M^){Pf+M^)' ^2.9) 



/ 



where the curly brackets denote symmetrization of the indices. Depending on the types of hadrons 
and currents involved, the resulting scattering amplitudes J^^ can describe a variety of different 
reactions such as Compton scattering, pion electroproduction, Nn or nn scattering, or crossed- 
channel processes such as pp annihilation into two photons or meson production. 

We have now outlined a systematic approach to compute various hadron properties from their 
underlying non-perturbative substructure in QCD. Its input is provided by the quark propagator S 
and the qq, qq and qqq kernels. Once these quantities are determined, no further model input is re- 



quired: we can selfconsistently solve the bound-state equation ( |2.3[ ) to obtain a hadron's wave func- 
tion and mass, solve the inhomogeneous BSE ( |2.7D for the vertex, and combine them to calculate 
hadron form factors and scattering amplitudes. The Dirac-Lorentz structure of the wave functions, 
vertices and current matrix elements is fully determined from Poincare covariance. By implement- 
ing the complete 'operator basis' in each case, its momentum-dependent Lorentz-invariant dressing 
functions (the 'form factors') are obtained from the equations described above. 

In order to proceed, we have to specify a truncation procedure. The apparent problem is the 
lack of information on the kernels which, according to the reasoning so far, encode the information 
from QCD's Green functions that is relevant for hadron physics. Omitting the term /Tp] in Eqs. ( ^!2| ) 
and yields the covariant Faddeev equation together with its corresponding currents [Q, |5|]. 
They trace the binding mechanism of three quarks in a baryon to its quark-quark correlations. 
Moreover, the simplest ansatz for the qq (and qq) kernel is the rainbow-ladder kernel: 

K = zl^-^^^Try^®y\ (2.10) 

where Tj^^ = 8^^ — k^k^ /k^ is a transverse projector with respect to the gluon momentum k. Its 
implementation in the Faddeev equation yields, by iteration, all dressed-gluon ladder exchanges 
between quark pairs. Implemented in the current, only the direct couplings to the quarks, i.e., the 
first two terms in Fig. ^ survive. 

The rainbow-ladder kernel satisfies vector and axialvector Ward-Takahashi identities which 
ensure electromagnetic current conservation and the Gell-Mann-Oakes-Renner and Goldberger- 
Treiman relations at the hadron level [|ll|, ^. Through these identities, Eq. ( 2.10 ) also determines 



the kernel of the quark DSE whose solution is numerically straightforward. The quark-gluon vertex 
is thereby reduced to its vector structure ~ y^, and its nonperturbative dressing, together with that 
of the gluon propagator, is absorbed in an effective interaction a{k^) which is the remaining model 



input. We employ the Maris-Tandy ansatz of Ref. [12] which reproduces the one-loop logarithmic 
running at large gluon momenta and features a Gaussian bump in the infrared that generates dy- 
namical chiral symmetry breaking. Different parametrizations have been recently tested and yield 



similar results for a range of hadron properties [13-15]. 
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Figure 3: Quark-mass dependence of nucleon static electromagnetic properties compared to lat- 
tice results. Left panel: squared isovector Dirac radius {r\)^. Right panel: isovector and isoscalar 
anomalous magnetic moments K^' and k" in units of nuclear magnetons. Stars denote the experi- 
mental values. Figure adapted from Ref. [S]. 



3. Nucleon electromagnetic form factors 

The nucleon's electromagnetic current is expressed by two dimensionless form factors: the 
Dirac and Pauli form factors Fi [Q^) and f2(G^). or the Sachs form factors Ge{Q^) and Gm{Q^) as 
their linear combinations: Ge = Fi — / (4M^)F2 and Gm = F1+F2. The current matrix element 
is given by 

Af, (3.1) 

where Q = Pf — Pi is the photon momentum and A^^^ — \ +^',/) positive-energy projectors. 
In the static limit one retrieves the nucleons' anomalous magnetic moments K = ^2(0) as well as 
their Dirac and Pauli radii r\ = — 6Fj'(0) and r\ = —6^2(0) 7^2(0). The isoscalar (isovector) form 
factors are the sum (difference) of proton and neutron form factors: F-'^ = Ff ±F". 

Results for the pion-mass dependence and 2^— evolution of various nucleon electromagnetic 
form factors are shown in Figs. |3| and ^ The bands correspond to a variation of the infrared prop- 
erties in the quark-gluon interaction a{k^) and measure the model uncertainty. As anticipated, the 
absence of pion-cloud contributions in the chiral and low-momentum region is recovered in the 
results. All form factors are in reasonable agreement with experimental data at larger momentum 
transfer where the nucleon is probed at small length scales and the pion cloud becomes irrele- 
vant. Missing structure mainly appears in the low-momentum region 2^ < 2 GeV^. The calculated 
charge radii, such as the isovector Dirac radius in the left panel of Fig. ^ underestimate their ex- 
perimental values but converge with lattice data at larger quark masses. Pion loops would increase 
the charge radii toward the chiral limit where they would diverge. 

Chiral effective field theory predicts that leading-order chiral corrections to proton and neutron 
anomalous magnetic moments carry an opposite sign; their magnitude is therefore enhanced in the 
isovector combination K^' = — k" and cancels in the isoscalar case k" = + k". The isoscalar 



J^" = iA 



f 



2M^ 
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Figure 4: Nucleon electromagnetic Sachs form factors as functions of the photon momentum trans- 
fer and in comparison with experimental data. The proton's electric form factor in the top left panel 
is normalized by the standard dipole. Figure adapted from Ref. [|]]. 



magnetic moment is quite accurately reproduced by the Faddeev calculation: jc* = —0.12(1), com- 
pared to the experimental value k^^^ = —0. 12 [||]. The calculated values of k:* and K^' correspond to 
an underestimation of 20% — 30% in the proton and neutron magnetic moments Gj^"(0), visible in 
the bottom panels of Fig. Q Another example is the neutron electric form factor G|(2^) in Fig. ^ 
which agrees with recent measurements at larger but misses the characteristic bump at low Q^. 
These observations suggest to identify the rainbow-ladder truncated nucleon with the 'quark core' 
in chiral effective field theories. 

The large-2^ behavior of form factors is of great theoretical and experimental interest as 
well. The experimental falloff of the proton's form factor ratio G^/G^ has been attributed to 
orbital angular-momentum correlations in the nucleon wave function which modify the perturbative 
scaling behavior and entail a zero crossing in G^(2^). Quark orbital angular momentum in terms 
of s, p and d waves appears in the Dirac-Lorentz structure of the nucleon's rest-frame Faddeev 
amplitude. While nucleon and A baryons are dominated by s waves, p waves play an important 
role as well: they contribute ~ 30% to the nucleon's canonical normalization and diminish only 
slowly with increasing current-quark masses. The contribution from d waves, on the other hand, is 
below 1%. At large Q^, the form-factor results from the Faddeev calculation become sensitive to 



7 



Baryon form factors from Dyson-Schwinger equations 



the numerics; nevertheless, a decrease of G| compared to the dipole form is visible in Fig. ^ and 
implies a zero crossing as well. 

Another remark concerns the timelike behavior of the form factors and the vector-meson dom- 
inance property which is a consequence of the underlying dynamics. The electromagnetic current 
is microscopically represented by the quark-photon vertex which can be separated in two terms: a 
Ball-Chiu part that satisfies electromagnetic gauge invariance, and another purely transverse term 
that includes vector-meson poles in the = 1— channel [|T§ |l7]|. Since the rainbow-ladder 
truncation does not dynamically develop hadronic decay widths, the poles that are generated in 
the self-consistent calculation of the quark-photon vertex are timelike and real. The decomposi- 
tion into 'Ball-Chiu' and 'p-meson' contributions can be made in all electromagnetic hadron form 
factors which therefore possess poles at = —nip and further 1 excited-state locations. The 
transverse term is negative at spacelike and, in the case of electric form factors, vanishes at 
= 0, i.e., the Ball-Chiu part alone satisfies charge conservation G^{0) = 1. The p— meson term 
contributes roughly ~ 50% to the nucleon's squared charge radii throughout the current-mass range 
but has only a minor impact on its magnetic moments whose overall contribution comes from the 
Ball-Chiu term. 

We note that a reduction of the Faddeev equation to a quark-diquark description, where scalar 
and axialvector diquark correlations are calculated from the same quark-gluon input, yields quite 
similar results for the form factors [18]. The model dependence is however larger, especially at 
large Q^, and the corresponding bands in Fig. ^ become sizeable; cf. also Fig. ^ below. Neverthe- 
less, these results imply that the interaction of quarks with scalar and axialvector diquarks provides 
the overwhelming contribution to the nucleon's binding. 



4. Nucleoli axial form factors 

In Ref. ||6|], the Faddeev approach has been further applied to compute the nucleon's axial and 
pseudoscalar form factors. The respective current matrix elements are specified by the axial form 
factor Ga{Q^), the induced pseudoscalar form factor Gp{Q^), and the pseudoscalar form factor 

/5=G5(e')A;/75A^ (4.1) 

Their microscopic decomposition in the Faddeev framework is identical to Fig. ^ except for the 
type of qq vertices that are involved: the structure that enters the self-consistent calculation 
of the quark-photon vertex is replaced by y^y^ and 75, respectively. Again, the pole structure of 
the resulting axial and pseudoscalar vertices allows to extract information on the timelike behav- 
ior and identify the relevant scales in the form factors. Ga is dominated by the 1++ axialvector 
meson ai(1260) and its excitations whereas Gp and G5 are governed by the pion pole. The pion- 
nucleon form factor Gj^nn is the residue of G5 at the pion pole and thus related to the 71 (1300) 
and further 0^+ excitations. The Goldberger-Treiman relation Ga(0) = /^G;rww(0)/M]v follows 
as a consequence of the axialvector Ward-Takahashi identity and analyticity which are satisfied 
microscopically. 



•'5 
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GA(e')7''+Gp(G') 
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Figure 5: Left panel: Quark-mass dependence of the nucleon's axial charge gA, compared to lattice 
results and the chiral expansion of Ref. [19]. Right panel: 2^— evolution of the axial form factor 



Ga(2^), compared to lattice data and the experimental dipole form. Figure adapted from Ref. 



The (isovector) axial and pseudoscalar form factor results exhibit various similarities with 
their electromagnetic counterparts, see Fig. ^. The axial charge gA = Ga{0) underestimates the 
experimental value by 20% — 25%; it falls below recent lattice data in the low quai^k-mass region 
and approaches the chiral expansion at lai^ger pion masses. On the other hand, Ga{Q^) is consistent 
with the phenomenological dipole form at larger Q^. Analogous results are obtained for the remain- 
ing pseudoscalar- form factors. This suggests once again that these features ai^e signals of missing 
pion-cloud effects. Such an interpretation was also proposed to explain the volume dependence of 
lattice results for gA [^]. 



5. Electromagnetic N ^ A transition 



Finally, the approach can be applied for the calculation of A(1232) and N ^ A transition form 
factors as well. Since a solution for the A bound-state amplitude from the Faddeev equation has be- 



come available only recently []21[], we will restrict our discussion to the quark-diquark model. The 
derivation that leads to the diagrams in Fig. ^ yields analogous expressions in the quark-diquark 
approach [22, ^], where the diquai^k ingredients can be computed self-consistently from the same 
quark-gluon input. Form-factor results in that framework exist for nucleon and A electromagnetic 

and in the following we will sum- 



form factors [ ^8| , p3| ] and the AA'^ 7i pseudoscalar transition [|2£ 
marize recent results for the electromagnetic A'^A 7 transition [^. 

The A'^Ay transition is characterized by the three Jones-Scadron form factors G]^(2^), 
and G^{Q^) which are related to the pion electroproduction multipole amplitudes at the A— resonance 
position [ ^ The respective cuiTcnt J^'P is decomposed as follows: 



-bpp/irs 



1(0 

2X7 



n 

CO 



GcQ"K^ 



A1 



(5.1) 



Instead of the incoming and outgoing momenta P,- and Pf, we used the orthonormal four-momenta 
and = Pj^ , where a hat denotes normalization, P = {Pf + Pi)/2 is the average momentum 
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Figure 6: 2^— dependence of the electric and Coulomb quadrupole form-factor ratios Rem and Rsm 
compared to experimental data. Figure adapted from Ref. 



and Pt its component transverse to the photon momentum. The Rarita-Schwinger projector for the 
A— baryon reads 

Pj« = A+ - 7^/ ) r/^" , (5.2) 



and the transverse projectors Tp^^ , Tq^ and T^^ are defined in the same way as in the gluon case 



below Eq. (2.10). The remaining dimensionless variables in (5.1) are given by 



2(M2+m2)' 2(M|+M2) ' • V + • V2V Mm 



The A'^Ay transition has been accurately measured over a wide momentum range [g6[ |27| ]. 
It is dominated by a magnetic dipole transition which, in a quark-model picture, amounts to the 
spinflip of a quark and is related to the form factor The remaining electric and Coulomb 

quadrupole form factors are much smaller and expressed by the ratios Rem{Q^) and Rsm{Q^) which 
encode the deformation in the transition. In non-relativistic quark models, non-zero values for these 
ratios would require J— wave components in the nucleon and A wave functions. On the other hand, 
the analysis of pion electroproduction data via dynamical reaction models suggests that Rem and 



Rsm are almost entirely dominated by the pion cloud [28]. 

In contrast, the quark-diquark results which are plotted in Fig. ^ reproduce the experimental 
data for Rem and Rsm quite well, even without the inclusion of pion-cloud corrections. In the case 
of Rem^ this behavior originates from wave contributions in the nucleon and A wave functions 
which are a consequence of Poincare covariance. The removal of p waves results in a ratio that 
is overall positive and grows with increasing Q^, cf. Fig. ^ with a trend towards the perturbative 
prediction Rem — ?> 1 for 2^ — )■ oo |27]. The impact of d waves is almost negligible. 



On the other hand, the result for the magnetic dipole transition form factor G\i{Q^) in Fig. ^ 
follows the characteristics of the previously discussed magnetic and axial form factors: it agrees 
with experimental data at larger and underestimates them by ~ 25% at = 0. This is consis- 
tent with the quark-model result and the expected behavior of the pion cloud from coupled-channel 
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Figure 7: Left panel: 2^— dependence of the magnetic dipole transition form factor G\f{Q^) com- 
pared to experiment. Right panel: decomposition of Rem{Q^) according to the orbital angular- 
momentum content in the nucleon and A wave functions. Figure adapted from Ref. 

analyses. Moreover, neither G^j nor Rsm are sensitive to the addition of p and d waves but domi- 
nated by 5— wave elements alone. 

6. Conclusions and outlook 

We have discussed several recent nucleon and A form factor results in the Dyson-Schwinger 
approach, obtained either directly from the covariant Faddeev equation or in a quark-diquark sim- 
plification. All calculations share the same quark-gluon input and the results display consistent 
features. Quark-quark correlations, which are mediated by a rainbow-ladder gluon-exchange inter- 
action, can account for the overall properties of the nucleon and A quark core and justify a quark- 
diquark picture for these baryons. Dynamical chiral symmetry breaking and Poincare covariance 
have important consequences for the behavior of the form factors. Their timelike structure is dom- 
inated by meson poles in the underlying quark-antiquark vertices. The admixture of quark orbital 
angular momentum via p waves, even in 5— wave dominated ground states such as the nucleon 
and A— baryon, is crucial for the NAy electric quadrupole form factor and the large-Q^ behavior 
of electromagnetic form factors. The main missing ingredients in a rainbow-ladder approach are 
pion-cloud contributions at low momenta and small pion masses. 

The combination of Dyson-Schwinger and covariant bound-state equations provides valuable 
tools for investigating the internal structure of hadrons. Its applications are still at an early stage, 
and it is desirable to extend the framework to study more sophisticated systems and reactions such 
as baryon excitations and nucleon-to-resonance transition form factors, virtual Compton scattering, 
pion electroproduction, pion-nucleon scattering, or timelike pp annihilation processes. At the same 
time, these efforts must be complemented by technical improvements, such as residue calculus to 
provide kinematic access to truly large Q^, or the implementation of pion-cloud corrections and 
hadronic decay channels via truncations beyond rainbow-ladder. 



11 



Baryon form factors from Dyson-Schwinger equations 



Acknowledgements 

I would like to thank C. S. Fischer and D. Nicmorus for useful discussions. This work was 
supported by the Austrian Science Fund FWF under Erwin-Schrodinger-Stipendium No. J3039, 
and the Helmholtz International Center for FAIR within the LOEWE program of the State of Hesse. 

References 

[1] R. Alkofer and L. von Smekal, Phys. Kept. 353 (2001) 281. 

[2] C. S. Fischer, J. Phys. G32 (2006) R253. 

[3] L. Chang, C. D. Roberts, and P. C. Tandy, Chin. J. Phys. 49 (201 1) 955. 

[4] G. Eichmann, R. Alkofer, A. Krassnigg, and D. Nicmorus, Phys. Rev. Lett. 104 (2010) 201601 . 

[5] G. Eichmann, Phys. Rev. D 84 (201 1) 014014 . 

[6] G. Eichmann and C. S. Fischer, arXiv:l 1 1 1.2614 [hep-ph]. 

[7] G. Eichmann and D. Nicmorus, arXiv; 1 1 12.2232 [hep-ph]. 

[8] A. N. Kvinikhidze and B. Blankleider, Phys. Rev C 60 (1999) 044004 . 

[9] M. Oettel, M. Pichowsky, and L. von Smekal, Eur Phys. J. A8 (2000) 251. 

[10] G. Eichmann and C. S. Fischer, ai-Xiv:l 1 1 1.0197 [hep-ph]. 

[11] P Maris, C. D. Roberts, and P C. Tandy, Phys. Lett. B 420 (1998) 267 . 

[12] P Maris and P C. Tandy, Phys. Rev C 60 (1999) 055214 . 

[13] M. Blank, A. Krassnigg, and A. Maas, Phys. Rev D 83 (201 1) 034020 . 

[14] S.-X. Qin, L. Chang, Y.-X. Liu, C. D. Roberts, and D. J. Wi\s.on,Phys. Rev C 84 (201 1) 042202 . 

[15] H. Sanchis-Alepuz, R. Alkofer, G. Eichmann, and R. WilHams, arXiv:l 1 12.3214 [hep-ph]. 

[16] J. S. Ball and T-W. Chiu, Phys. Rev D 22 (1980) 2542 . 

[17] P Maris and P C. Tandy, Phys. Rev C 61 (2000) 045202 . 

[18] G. Eichmann, PhD Thesis, University of Graz, arXiv:0909.0703 [hep-ph]. 

[19] M. Procura, B. U. Musch, T R. Hemmert, and W. Weise, Phys. Rev D 75 (2007) 014503 . 

[20] S. Ohta, AIP Conf Proc. 1354 (201 1) 148. 

[21] H. Sanchis-Alepuz, G. Eichmann, S. Villalba-Chavez, R. Alkofer, Phys. Rev D 84 (201 1) 096003 . 

[22] M. Oettel, G. Hellstern, R. Alkofer, and H. Reinhardt, Phys. Rev C 58 (1998) 2459 . 

[23] D. Nicmorus, G. Eichmann, and R. Alkofer, Phys. Rev D 82 (2010) 1 14017 . 

[24] V. Mader, G. Eichmann, M. Blank, and A. Krassnigg, Phys. Rev D 84 (201 1) 034012 . 

[25] H. F Jones and M. D. Scadron, Ann. Phys. 81 (1973) 1-14 . 

[26] V. Pascalutsa, M. Vanderhaeghen, and S. N. Yang, Phys. Rept. 437 (2007) 125. 

[27] I. G. Aznauryan and V. D. Burkert, Prog. Part. Nucl. Phys. 67 (2012) 1-54. 

[28] B. Julia-Diaz, T.-S. H. Lee, T Sato, and L. C. Smith, Phys. Rev C 75 (2007) 015205 . 



12 



